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ABSTRACT 

The  axi symmetric  vortex  motion  of  a  viscous  fluid  ever  a  flat 
surface  is  investigated.  Paradoxical  results  obtained  by  Moore  and 
Goldshtik  are  resolved.  Guided  by  a  useful  extension  or  prandtl's 
boundary  layer  theory  an  appropriate  similarity  transformation  is 
found  which  reduces  the  Navier-Stokes  equations  to  three  ordinary 
differential  equations.  The  solution  of  the  remaining  boundary  value 
problem  reveals  significant  properties  of  vortex  flows  and  explains 
several  phenomena  observed  in  hurricanes.  A  complete  numerical 
result  is  displayed. 
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1 .  Introduction 

In  recent  years  several  investigations  have  been  devoted  to  the 
mathematical  description  of  the  steady  vortex  motion  of  a  viscous  fluid 
which  is  governed  by  the  Navier-Stokes  equations.  Taylor  [9]  ,  Cooke 
[2],  Goldshtik  [3 1 ,  and  Long  [4]  investigated  the  vortex  motion  normal 
to  a  flat  surface,  which  is  asymptotic  to  a  coaxial  potential  vortex 
at  large  distances  from  the  surface  and  the  core  of  the  vortex.  Taylor 
and  Cooke  excluded  the  core  of  the  vortex  and  presented  certain  solu¬ 
tions  of  the  corresponding  classical  boundary  layer  equations  by  apply¬ 
ing  the  method  of  Kdrman-Pohlhausen .  Moore  [5]  showed  that  the 
axisymmetric  vortex  flow  considered  is  nonasymptotic  to  certain  so- 
called  "similarity  solutions"  of  the  boundary  layer  equations  used  by 
Taylor  and  Cooke.  Nevertheless,  Long  ignored  the  presence  of  the  sur¬ 
face  and  constructed  a  certain  similarity  solution  of  the  classical 
boundary  layer  equations ,  which  presumably  describe  the  vortex  flow 
approximately  in  the  core.  Goldshtik  reexamined  the  problem  and 
attempted  to  integrate  the  exact  Navier-Stokes  equations  by  a  similarity 
solution  of  Moore's  type.  This  interesting  investigation  led  to  the 
paradoxical  results  presented  in  [3]  , 

Since  the  paradoxical  results  of  Moore  and  Goldshtik  appear  to 
be  physically  questionable,  the  same  vortex  flow  problem  will  be  rein¬ 
vestigated  in  the  present  paper.  It  will  be  shown  that  the  dis¬ 
crepancies  discovered  by  Moore  and  Goldshtik  are  due  to  inconsistent 
and  incomplete  boundary  data  used  to  determine  a  unique  solution  of 
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the  differential  equations  considered .  Moore's  nonexistence  result 
is  also  a  consequence  of  Prandtl's  boundary  layer  theory  which 
destroys  the  elliptic  character  of  the  Navier-Stokes  equations  without 
examining  the  existence  of  a  solution  to  the  remaining  parabolic 
boundary  value  problem. 

On  the  basis  of  a  specified  physical  problem ,  complete  and 
consistent  boundary  data  will  be  found  in  order  to  define  a  useful 
solution  of  the  Navier-Stokes  equations.  Guided  by  an  extended 
boundary  layer  theory ,  which  preserves  the  full  elliptic  Navier- 
Stokes  equations ,  a  set  of  ordinary  differential  equations  will  be 
derived  by  an  appropriate  similarity  transformation.  This  trans¬ 
format  ion  is  based  upon  the  concept  of  the  limiting  line  of  a 
boundary  layer  the  existence  of  which  is  suggested  by  the  extended 
boundary  layer  theory. 

The  investigations  are  supported  by  a  complete  numerical  solu¬ 
tion,  which  has  been  established  by  means  of  the  Runge-Kutta  method . 
All  results  appear  to  be  physically  plausible  and  permit  the 
explanation  of  several  phenomena  observed  in  hurricanes  outside 
their  cores. 
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2 .  Definition  of  a  Solution  of  the  Navier-Stokes  Equations 

In  the  cylindrical  coordinate  system  (r ,  tp,  z)  let  (u,  v,  w) 
denote  the  corresponding  velocity  vector  of  an  axisymmetric  vortex 
flow  over  the  solid  surface  z  «*  0.  If  V,  p,  and  p  are  the  constant 
kinematic  viscosity,  the  constant  density,  and  the  variable  pressure 
of  the  fluid,  the  flow  is  governed  by  the  Navier-Stokes  equations 


uur  +  wu2 


uv..  +  w v„  +  — 
r  z  ^ 


uw  +  ww 
r  z 


f  |  »r  *  ♦  (f)t  ♦  ««] 

+  ”[V«  +  (f)r  "•''«] 

‘ipz+  «frt+  K  +  “«J 


r  +  (^>2  "  0 


(1) 


(2) 


(3) 

(4) 


Since  the  nonlinear  elliptic  partial  differential  equations  (1) 
through  (4)  are  singular  at  the  axis  r  ■  0,  it  is  not  possible  to 
specify  a  solution  in  the  region  (r  s  0,  z  5  0)  by  an  arbitrary  set 
of  boundary  data  and  singularity  conditions .  In  order  to  find  a 
consistent  and  complete  set  of  boundary  values ,  it  is  helpful  to 
examine  the  characteristic  properties  of  the  following  physical 
flow  model  (see  Figure  1) . 

The  vortex  flow  over  the  surface  at  z  ■  0  may  be  produced 
by  a  very  long  rod  of  small  diameter  d  which  is  rotating  around  its 
axis  at  r  •  0  with  the  angular  velocity  w.  If  all  friction  forces 
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at  the  surface  at  z  =  0  were  eliminated,  then  the  rod  would  produce  the 

potential  vortex  flow 

u  *  0»  v  ml  »  w  £  0,  l  £  *  (5) 

r  p  r3 

with  the  vortex  strength 

r  -  »(f)a  ,  <«> 

where  an  additive  pressure  Constant  is  neglected.  In  order  to  investi¬ 
gate  exclusively  the  influence  of  the  friction  forces  at  the  surface  z  =  0 
on  the  potential  vortex  (5),  it  is  necessary  to  prevent  any  additional 
exterior  disturbances  Of  the  flow.  This  condition  may  be  fulfilled  by 
providing  the  rod  with  a  flexible  surface  which  glides  freely  with  the 
flow  in  the  axial  direction.  For  the  rigorous  mathematical  model  the  rod 
must  extend  from  z  =  0  to  z  =  «.  The  diameter  Of  the  rod  must  Shrink 
to  zero,  while  the  angular  velocity  increases  to  infinity  so  that  the 
vortex  strength  (6)  remains  a  constant.  Since  the  fluid  is  to  be  at 
rest  at  large  distances  from  the  rod,  adequate  sources  and  sinks  of 
equal  strength  must  be  located  at  large  distances  from  the  surface. 

A  circular  sink  around  the  axis  r  =  U  is  necessary  at  z  =  »  to  absorb 
the  fluid  which  is  ejected  by  the  vortex  motion  from  the  boundary  layer 
at  the  surface,  This  sink  must  change  at  some  r  =  r0  to  an  annular 
source  of  equal  strength,  which  supplies  the  boundary  layer  with  the 
necessary  fluid  in  return. 

With  this  description  of  the  physical  flow  model,  it  is  now 
easy  to  determine  consistent  singularities  and  complete  boundary  data 
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which  define  a  useful  solution  of  the  Navier-Stokes  equations.  At  the 
surface  and  at  large  distances  from  the  rod  the  solution  is  subject  to 
the  "tegular  boundary  conditions" : 

*  >  q  }  :  u  »  0,  v  -  0,  w  «  0  (7) 


r  •  “  h  u 

z  <  »  J 


0, 


0, 


(8) 


At  the  rotating  rod  and  at  large  distances  from  the  surface  the  solu¬ 


tion  is  determined  by  the  following  "singular  boundary  conditions" : 


r  -  9  }  :  u  -  0  , 

i  >  0  J 


•  1 


TA  log  y0- 


(9) 


r  <  •  c  :  u  -»  0  ,  —  1  ,  -  l  ,  (10) 

»  ~  •  -  r  rA  log  fa 

where  T  and  A  are  constant  parameters  at  one's  disposal.  The  constant 

tp  must  be  determined  simultaneously  with  the  solution.  Mo  conditions 
are  Imposed  at  the  two  singular  points  (r  «  0,  z  =  0)  and  (r  ■  »,  z  *  ») , 
The  consistency  of  the  foregoing  singularities  (9)  and  (10)  may 
be  demonstrated  on  the  basis  of  the  following  principles : 

(i)  The  singularities  required  for  the  primary  tangential  flow  v  are 
admissible,  because  the  Mavier-Stokes  equations  (1)  through  (4) 
yield  the  potential  solution  (5)  which  has  the  same  singularities. 

(it  The  sink  and  source  distribution  admitted  for  the  secondary  axial 
flow  w  is  compatible  with  the  singularities  of  the  primary  tan¬ 
gential  flow,  as  the  Navier-Stokes  equations  yield  the  potential 
solution 

F 

u  «  0  ,  v*®,  w  *  TA  log  (11) 
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with  the  same  singularities .  Thus  the  exchange  of  fluid  between 
the  boundary  layer  at  z  =  0  and  the  Sink  and  source  at  z  *  »  Can 
proceed  freely  without  friction  forces  acting,  since  the  undis¬ 
turbed  flow  (5)  is  not  regular  at  r  =  ®,  the  amount  of  fluid 
revolved  by  the  vortex  flow  depends  upon  the  radius  r  =  ®  of  the 
surface  at  z  =  0.  This  indicates  the  freedom  of  the  parameter  A 
which,  vice  versa,  determines  the  radius  of  the  surface,  i.  e. ,  it 
defines  the  exact  meaning  of  the  limit  r  =  ®  (see  [8]) . 

In  this  connection  it  is  significant  to  compare  the  boundary 
data  posed  by  M.  A.  Goldshtik  for  the  same  physical  problem  (see  [3] 
and  also  [2,  5,  9 ]) ,  Goldshtik  Subjected  the  flow  to  the  same  regular' 
boundary  values  (7)  and  (8) ,  to  Which  he  added  the  singular  boundary 
conditions 


r 

z 


hi 

i* 


i-i 

p 


(12) 


r  -  0  :  . 
z  >0  J  * 


u  -  0  , 


w  -•  w0(z)  (bounded  1) 


(13) 


Since  these  conditions  were  incomplete,  he  required  later  a  special 
"similarity  property"  for  the  solution,  by  which  he  tacitly  added  the 
missing  data 


r  >  Q 

Z  -  05 


W  -  0  ,  (»  -  0)  , 


(14) 


r  -  0  i  .  rv 
z  >  0  J  '  r 


1  , 


(15) 
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Without  going  into  any  process  of  integration,  it  is  easy  to 
show  the  inconsistency  of  these  conditions  Which  led  Goldshtik  to  hit 
"paradoxical  results".  His  special  similarity  assumption  requires  the 
vortex  flow  to  produce  a  sort  of  source  in  the  axial  flow  w  at  the 
point  (r  ■  0,  z  =  0) .  This  source  is  characterized  by  an  unbounded 
axial  velocity  w  which  must  remain  unbounded  along  the  entire  axis 
r  =  0,  where  the  tangential  and  radial  flows  are  at  their  asymptotic 
frictionless  state.  In  fact,  Goldshtik' s  existence  proof  for  small 
Reynolds  numbers  is  based  on  the  second  condition  (13),  which  is  not 
fulfilled  by  his  Solution  (see  [3,  page  927  ^)? 

w*  ..c-.  .  (16) 

*  — v  z 

It  may  be  mentioned  that  the  authors1  attempt  to  Solve  Goldshtik' s 
equations  (1.3)  through  (1.8)  tinder  the  boundary  conditions  (1.7) 
and  (1.8)  (see  [3])  failed  for  Small  and  large  Reynolds  numbers 
because  of  nonconvergence  of  the  Runge-Kutta  method.  The  numerical 
results  demonstrated  very  clearly  that  the  friction  forces  at  the 
Surface  are  not  sufficient  to  produce  a  source  of  the  strength  (lb) 
in  the  axial  flow.  Further  exterior  forces  are  needed,  to  reduce  the 
axial  velocity  from  w  *  »  to  w  e  0  along  the  axis  r  =  0,  This  obser¬ 
vation  led  to  the  introduction  of  the  logarithmic  singularity  (11) 
which  is  free  of  effective  friction  forces. 
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3.  A  Reduction  of  the  Navier-Stokes  Equations 

For  an  approximate  integration  of  the  Navier-Stokes  equations  (1) 
through  (4)  under  the  boundary  conditions  (7)  through  (10) ,  it  is  useful 
to  simplify  the  singular  boundary  conditions  (9)  and  (10)  by  introducing 
dimensionless  quantities  by  the  following  conventions : 


u  ■ 


w  ■  TA  log  |n  w> 


(17) 


£  - 

P 


(18) 


where  R  may  be  called  the  Reynolds  number  of  the  flow.  After  carrying 
out  this  transformation  one  arrives  at  the  new  Navier-Stokes  equations: 

+  at  Ug  |»  -  iflr1  +  v3)  -  pt  p  +  *  #„)  -  vt:  (W) 

uvr  +  Ar  log  £a  WVZ  ‘  -  I[r(Vrr  +  Vzz)  -  Vr]  (20) 

r  R 

r  log  |a[uwr  +  Ar  log  WWZ]  -  UW 

Ur  +  Ar  log  Wz  ■ 
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The  transformed  boundary  conditions  are: 


r  >  01  j  o  *  0,  V  -  0,  W  «  0,  (23) 

z  «  0  J 

l  “  *}'  :  W  m  0,  V  -  0,  W  -  0,  (24) 

S  >2}  :  U  “  °*  v  “  l»  W  *  X*  (25) 

l*  j  i  «  -  °.  V  -  1,  W  -  1  .  (26) 


The  solution  of  this  equivalent  elliptic  boundary  value  problem  may  be 
guided  by  the  essential  hypothesis  of  Prandtl's  boundary  layer  theory  in 
the  following  form. 

Weak  Boundary  Layer  Assumption:  Friction  forces ,  which  are 

caused  by  the  nonslip  condition  at  a  solid  surface  in  a  primary  real 
flow  without  effective  friction  forces,  are  essentially  acting  only 
within  a  "boundary  layer"  that  is  bounded  by  a  certain  "limiting  line." 
The  asymptotic  flow  beyond  the  boundary  layer  is  the  primary  flow  on 
which  may  be  superposed  a  secondary  flow  without  effective  friction 
forces . 

It  is  significant  to  note  that  the  "extended  boundary  layer 
theory"  based  on  the  weak  assumption  preserves  the  character  of  the 
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Navier-Stokes  equations .  Thus  the  validity  of  existence  and  unique¬ 
ness  theorems ,  which  are  known  for  elliptic  boundary  value  problems, 
remains  untouched  by  this  theory .  This  most  important  requirement  of 
admissible  simplifications  of  mathematical  problems  is ,  however,  ignored 
by  the  complete  "classical  boundary  layer  theory."  Indeed ,  the  well- 
known  boundary  layer  assumptions,  which  neglect  certain  first  and 
second  order  partial  derivatives,  lead  to  parabolic  differential 
equations ,  for  which  the  existence  and  uniqueness  theorems  of  elliptic 
boundary  value  problems  do  not  hold  in  the  same  generality.  This  is 
reason  enough  to  abandon  such  risky  assumptions.  In  this  connection 
it  may  be  interesting  to  compare  other  simplifications  of  this  sort 
which  have  been  suggested  in  the  past  and  which  have  been  criticized 
in  a  similar  manner  (see,  for  instance, [l]) . 

The  additional  assumptions  of  the  classical  boundary  layer 
theory  may  also  be  abandoned  as  they  do  not  simplify  the  constructive 
solution  of  the  remaining  mathematical  problem  in  general.  This 
statement  will  be  verified  for  the  present  problem  by  the  following 
derivations  and  for  further  examples  by  other  investigations  (see  [8]), 
It  may  be  worthwhile  to  mention  that  the  extended  boundary  layer 
theory  introduced  above  leads  to  exactly  the  same  asymptotic  solution 
of  the  classical  flat  plate  problem  as  the  common  boundary  layer 
theory  without  any  mathematical  complications .  This  reveals  a  clear 
superiority  of  the  extended  boundary  layer  theory  over  the  classical 
theory  which  may  be  well  founded  in  special  cases. 
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In  order  to  find  an  approximate  solution  of  the  boundary  value 
problem  considered,  it  is  helpful  to  utilize  the  limiting  line  z  »  6(r) 
of  the  boundary  layer  the  existence  of  which  is  suggested  by  the 
extended  boundary  layer  theory .  For  an  accuracy  parameter  e,  which 
is  at  one's  disposal,  the  limiting  line  of  the  boundary  layer  along 
the  surface  at  z  =  0  may  be  characterized  by  the  following  flow  data: 

z  ■  6(r)  :  V  “  1  -  c  (27) 

*  i  4(r)  :U«0  ,  V  1  ,  W^l  .  (28) 

The  existence  of  a  line  z  =  6(r) ,  which  satisfies  all  conditions  (27) 
and  (28)  is  plausible  from  the  boundary  data  (23)  through  (26) . 

Since  V  is  a  solution  of  an  elliptic  differential  equation,  6(r)  must 
be  an  analytic  function  of  r  which  connects  the  two  singular  boundary 
points  (r  =  0,  z  **  0)  and  (r  ■  «#,  z  *  »)  for  any  accuracy  parameter 
0  <  «  <  1  (see  Figure  1) . 

According  to  the  conditions  (23)  through  (28)  it  is  useful  to 
transform  the  Navier-Stokes  equations  (19)  through  (22)  by  the 
"similarity  relation" 


which  maps  the  boundary  layer  0  ^  z  £  6(r)  onto  the  parallel  strip 
0  s  (  s  l,  The  boundary  lines  (r  >  0,  z  ■  0)  and  (r  <  ■»,  z  “  ®)  are 
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mapped  onto  the  corresponding  lines  (r  >  0,  £  ■  0)  and  (r  <  ®,  £  =  ®) , 

While  the  lines  (r  =  ®,  z  <  »)  and  (r  »  0,  z  >  0)  reduce  to  the  points 

(r  =  ®,  £  =  0)  and  (r  =  0,  £  =  ®) ,  the  singular  points  (r  =  0,  z  =  0) 

and  (r  =  ®,  z  »  ®)  stretch  into  the  lines  (r  *  0,  £  <  ®)  and  (r  «  ®,  £  >  0) . 

Thus ,  in  the  (f,  £)  plane  the  solution  Conr.r'.ierecl  is  not  ^overr.ed.  by 

any  boundary  data  at  r  =  0  and  r  =  ®  . 

After  evaluating  the  transformed  conditions  (23)  through  (28)  one  is 

tempted  to  seek  a  solution  of  the  corresponding  Navier-Stokes  equations 
which  is  independent  of  r.  Indeed ,  if  a  boundary  layer  exists ,  which 
satisfies  the  conditions  (27)  and  (28)  ,  then  all  partial  derivatives 
with  respect  to  r  vanish  at  £  *  0  and  are  almost  zero  at  £  =  1  at 
least  for  some  values  of  r,  for  instance,  around  r  »  r6  .  The  same 
remains  true  within  the  boundary  layer  0  £  £  §  1,  because  the  lines 
£  =  const .  may  be  considered  as  first  order  approximat ions  of  the 
lines  of  constant  velocity  (U,  V,  W)  around  r  *  r 0  .  Thus ,  the  exact 
solution  of  the  problem  considered  must  yield  a  first  order  approxima¬ 
tion  which  is  independent  of  r  around  r  *»  r0  within  the  strip 

0  £  £  s  1. 

These  results  reveal  the  limits  of  the  classical  boundary  layer 
theory .  After  differentiating  the  essential  boundary  layer  assumption  (27) 


along  the  line  z.  ■  6(r)  one  arrives  at 


Thus, 


dr  dr  dz  dr 


(30) 


M  and  M>  display  the  same  asymptotic  behavior ,  if  |67(r)  |  has 
dr  dz 


a  positive  lower  bound.  Unless  r  can  be  confined  to  a  proper  vicinity 


NWL  REPORT  NO.  1835 

of  a  point  r  =  r0  at  which  6  /(r0)  *  0,  6  '(r)  may  even  grow  beyond 

any  upper  bounds,  if  the  accuracy  of  condition  (27)  is  increased. 

After  differentiation  of  equation  (30)  along  the  line  z  =  6(r) 
one  sees  in  the  same  manner  that,  in  the  Navier-Stokes  equations  (19) 
through  (22) f  partial  derivatives  with  respect  to  r  cannot  be  neglected 
up  to  second  order ,  unless  the  derivatives  6/(r)  and  6*(r)  vanish  at 
some  point  r  =  r0  .  Hence ,  the  classical  boundary  layer  theory  is 
consistent  for  an  "almost  parallel"  boundary  layer  which  has  a  limit¬ 
ing  line  that  assumes  the  slope  and  the  curvature  of  the  corresponding 
surface  in  the  region  of  interest.  This  condition  is  fulfilled  in 
the  classical  flat  plate  problem,  but  it  is  violated  in  the  present 
problem  as  the  following  derivations  will  show  (see  also  [8[), 

After  carrying  out  the  similarity  transformation  (25)  one 
arrives  at  the  reduced  Navier-Stokes  equations : 

(1  +  6/3Ca)U  +  [£(26 /2  -  66  '  +  i^)  +  ^(6'CU  -  Ar  log  |a  wftj 

+  R&u2  +  V  -  2P  -  CP]  -  0  (31) 

ra  6 

(1  +  6/2Cs)V  +  [£(26 ^  -  66"  +  M»)  +  M(6'£U  -  Ar  log  W) ]v  -  0  (32) 

w  +  XI  ^  •  5F  108  r^(l  +  *  C(26 ,p’  -  66  -  -p> 

+  |  log  |«[6  'CU  -  Ar  log  W]W  -  W  (33) 

•  i 

6  '£U  -  Ar  log  |ft  W  »  0  (34) 
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In  these  equations  the  partial  derivatives  of  U,  V,  W,  and  P  with 
respect  to  r  are  omitted  because  only  ^-dependent  solutions  are  of 


interest .  The  partial  derivatives  with  respect  to  Q  are  indicated 
by  a  dot ,  while  the  primes  denote  the  derivatives  of  the  limiting 
line  of  the  boundary  layer . 

The  continuity  equation  (34)  displays  very  clearly  that  the 
search  for  an  approximate  solution  of  the  boundary  value  problem 
considered  can  be  successful  in  the  present  form,  only  if 
6 '  (r)  yields  the  expansion 


A  Bro^j-j  log  +  * » •  (B  ^  0) 


(3$) 


for  instance  around  r  «  r0,  where  terras  of  second  and  higher  order  may 
be  neglected .  The  limiting  line  of  the  boundary  layer  has  then  the 


expansion 

6  -  B^)  £2  log  +  1 J  +  . . .  (B  4  0)  P6) 

around  r  =  r0  . 


Tlie  general  features  of  the  limiting  line  of  the  boundary  layer 
determined  by  equation  (36)  are  displayed  in  Figure  L.  It  is  remarkable 
to  note  that  the  boundary  layer  thickness  reaches  a  relative  maximum 
at  the  point  r  =  r0  where  the  axial  velocity  w  changes  its  direction. 
Hence  the  boundary  layer  thickness  must  assume  also  a  relative  minimum 
before  it  increases  rapidly  beyond  any  bounds .  This  surprising 
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phenomenon  is  physically  plausible,  as  the  positive  and  negative 
secondary  axial  flows  deform  the  boundary  layer  which  is  produced  by 
the  friction  forces  along  the  surface  at  z  ®  0. 

After  substituting  the  approximation  (36)  in  the  equations  (31) 
through  (34)  and  neglecting  terms  of  higher  order  one  obtains  the  fol¬ 
lowing  system  of  ordinary  differential  equations: 


BOi  -  "  *  0  •  (40) 

These  equations  may  be  simplified  by  the  introduction  of  the  stream 
function  G(Q  and  the  characteristic  number  a  which  are  defined  by 
the  relations 

U  =  -  G  ,  W  *  ?(G  -  CQ)  ,  c  ■  M£  .  (41) 

a  2 

This  substitution  leads  to  the  following  solution  of  the  problem 
considered : 

The  functions  G(Q,  V(»)f  and  P(C)  represent  an  approximate 
solution  of  the  Navier-Stok.es  equations  (1)  through  (4)  under  the 
boundary  conditions  (12)  through  (15) ,  which  is  valid  within  the 


boundary  layer  0  §  £  S  1  in  the  vicinity  of  r  «  r0 ,  provided: 
(A)  6(0  >  V(0  »  and  f(£)  fulfill  the  differential  equations 


G  +  0s  CG  -  (G8  +  v3  -  2P) 

4 

(42) 

V  +  0*CV  «  o 

(43) 

G(G  *  CG)  -  ^  P 

(44) 

Gf 


under  the  boundary  conditions 

£  -  0:  6  =  6,  G-0,  V  -  0  (45) 

£  «■  •:  G  »  G^,  V  -  l  ,  p  "  2  *  (A6) 

(B)  the  characteristic  number  a  satisfies  the  accuracy  condition 

C  ^  1*  V  -  1  -  «  ,  (47) 

(C)  G(0  end  P(0  satisfy  the  conditions 


C  -  1: 


G  0  ,  G  «0  , 
P  «  0 


(48) 


with  sufficient  accuracy, 

(P)  the  characteristic  constants  B  and  r0  of  the  limiting  line  of  the 
boundary  layer 


f 
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are  determined  by  the  relations 
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B  - 


2o 

ro  "  TT 


(50) 


This  solution  is  based  upon  the  "joining  properties"  (B)  and  (C) , 
which  constitute  the  mathematical  justification  Of  the  entire  boundary 
layer  theory.  The  exact  integral  (11)  of  the  Navier-Stokes  equations 
(1)  through  (4)  is  utilized  as  an  approximate  solution  outside  the 
boundary  layer  0  s  £  *  1  of  the  problem  under  consideration.  An 
approximate  solution  of  the  Navier-Stokea  equations  is  constructed 
within  the  boundary  layer  Which  yields  the  correct  boundary  data  at 
£  *  0.  Both  solutions  are  asymptotically  equal  for  large  values  of  £  . 
They  are  joined  along  the  limiting  line  of  the  boundary  layer  £  »  1 
under  a  small  violation  of  the  analyticity  of  the  exact  solution. 

This  violation  is  controlled  by  the  properties  (B)  and  (C) ,  which  is 
extended  to  all  derivatives  of  the  dependent  variables  that  occur  in 


the  Mavier •Stokes  equations. 

The  existence  of  an  integral,  which  fulfills  the  differential 
equations  (42) ,  (43) ,  and  (44)  and  the  boundary  conditions  (45)  and  (46) 
indicates  the  simplified  system  of  differential  equations 


•••  ••  R(Jr 

G  +  o3  £G  »  ~  (V3  -  1) 

V  +  o®CV  -  0  ,  pel 


(51) 

(52) 
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In  fast ,  (the  system  (51)  and  (52)  is  integrable  by  quadratures  and 
yields  a  solution  under  the  unaltered  boundary  conditions  (45)  and 
(46) .  Although  this  system  represents  no  approximate  substitute  lor 
the  system  (42),  (43) ,  and  (44) ,  it  retains  all  its  characteristic 
properties  at  both  ends  of  the  interval  0  s  Q  i  »  ,  Thus,  the  solution 
of  the  simplified  system  (51)  and  (52)  may  be  used  as  an  initial  inte¬ 
gral  for  an  appropriate  iteration  procedure  (see  section  4)  which  con¬ 
verges  to  the  solution  of  the  complete  system. 

The  existence  of  the  integral  defined  by  the  conditions  (A) 
through  (D)  justifies  the  general  boundary  layer  assumption  which 
guided  the  foregoing  derivations .  Furthermore,  it  demonstrates  the 
consequences  of  ignoring  the  restrictive  limits  of  Prandtl' s  boundary 
layer  theory.  If  all  assumptions  of  the  classical  boundary  layer 
theory  were  applied  then  the  differential  equations  (42) ,  (43),  and 
(44)  would  reduce  to  the  equations 


G  =  (62  +  v2  -  2P) 

(53) 

V  =  0 

(54) 

0*G(G  -  ;G)  “  P  , 

(55) 

which  have  no  solution  under  the  boundary  conditions  (45)  and  (46) , 

It  is  obviously  false  to  neglect  the  term  c*" ,G  against  all  other  terms 
of  the  differential  equations  (42)  and  (43)  no  matter  how  large  the 
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Reynolds  number  R  may  be  assumed .  It  is  exactly  this  term  which 
characterizes  the  boundary  layer  directly  at  the  surface  and  at  large 
distances  from  the  surface,  if  the  term  is  deleted,  then  the 
differential  equations  assume  a  totally  different  character  within 
and  outside  the  boundary  layer  and  no  longer  permit  a  proper  solution. 
As  was  pointed  out  above,  this  is  not  surprising  as  the  classical 
boundary  layer  theory  destroys  the  elliptic  character  of  the  Navier- 
Stofces  equations  without  examining  the  existence  of  a  proper  solution 
to  the  remaining  parabolic  equations. 
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4.  Characteristic  Properties  of  Vortex  Flows 

In  section  3  the  vortex  flow  problem  defined  in  section  2  has  been 
reduced  to  an  ordinary  boundary  value  problem  the  integral  of  which 
represents  an  approximate  solution  of  the  exact  problem.  Without 
complete  integration  of  the  remaining  boundary  value  problem  the 
foregoing  derivations  and  the  equations  obtained  display  significant 
phenomena  of  the  vortex  flows  considered . 

The  derivations  in  the  sections  2  and  3  led  to  two  important 
properties  of  vortex  flows ,  which  may  be  summarized  as  follows : 

Phenomenon  1 :  The  boundary  layer  along  a  flat  surface  normal  to 
an  axisymmetric  vortex  flow  may  be  considered  as  a  logarithmic  source 
and  sink  pair  for  the  secondary  axial  flow  outside  the  boundary  layer . 
The  radial  extension  r0  of  the  circular  axial  source  is  a  characteristic 

flow  parameter  in  addition  to  the  vortex  Reynolds  number  R  ■  t/v  , 
Phenomenon  2;  The  secondary  axial  flow  deforms  the  monotonic 
structure  of  the  friction  boundary  layer  along  the  flat  surface  as 
sketched  in  figure  1 . 

The  equation  (43)  ,  which  is  integrable  by  the  error  function  leads 
to  the  following  properties  of  the  vortex  flow  considered. 

Phenomenon  3 :  Within  the  accuracy  of  a  first  order  approximation 

the  dimensionless  velocity  (U,  V,  W)  is  a  vector  function  of  the 
dimensionless  variable  £  =  z/6(r),  such  that: 
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(a)  the  tangential  velocity  V(0  is  independent  of  both  characteristic 
flow  parameters  R  and  r0  , 

(b)  the  radial  and  axial  velocities  u(Q  and  w(Q  depend  solely  upon 
the  Reynolds  number  R  of  the  vortex  flow. 

Phenomenon  4 :  Within  a  first  order  approximation  the  characteristic 
number  a  (a  «  2.5  for  a  relative  accuracy  er  •  1%)  is  independent  of 
both  characteristic  flow  parameters  R  and  r0 .  Hence  the  boundary  layer 
thickness 


z  -  fi(r)  «  a  ^  ^  j  ^1+  2  log  (56) 

is  independent  of  the  Reynolds  number  R  of  the  vortex  flow. 

The  statement  (a)  of  the  phenomenon  (3)  is  an  obvious  consequence 
of  the  integral 

yr 

V(0  -  erf  (-=£)«-—  f  e-t3dt  ,  (57) 

w  i 

which  solves  equation  (43)  under  the  corresponding  boundary  condi* 
tions  (45)  and  (46) .  The  equations  (47)  and  (57)  lead  directly  to 
a  constant  number  <r  as  was  stated  in  phenomenon  4.  From  this 
result  one  deduces  easily  the  statement  (b)  of  phenomenon  3  by 
examining  the  equation  (42)  together  with  the  corresponding  boundary 
conditions  (45)  and  (46) , 
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In  this  connection  a  significant  analogy  between  the  present 
problem  and  the  classical  flat  plate  problem  may  be  pointed  out. 

In  both  cases  the  decay  of  the  corresponding  boundary  layers  along 
flat  surfaces  is  of  the  same  strength,  i.  e. ,  the  disturbances  of  the 
primary  velocities  decay  essentially  with  the  same  logarithmic  order 
two.  This  analogy  Seems  to  indicate  a  very  general  property  of 
boundary  layers  along  solid  surfaces  (see  also  [8 ~).  This  is 
physically  plausible  as  the  occurrence  of  a  friction  boundary  layer 
along  a  solid  surface  represents  a  transport  phenomenon. 

For  an  investigation  of  further  phenomena  of  vortex  f lows ,  it  is 
helpful  to  examine  complete  solutions  of  the  remaining  boundary 
value  problems  defined  in  the  previous  chapter.  Such  explicit  solu¬ 
tions  may  be  obtained  through  numerical  integrations  after  introduc¬ 
ing  the  following  new  variables 

^  ,  g  *  SG  ,  and  h  =  2P  ,  (58) 

This  substitution  leads  to  the  transformed  system  of  differential 
equations 


"g  +  Te  =  f(gs  +  V-  -  h) 
V  ♦  TV  •  0 


g(g  *  T|g)  =  ^h 


(59) 

(60) 


(61) 
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which  must  be  integrated  under  the  boundary  conditions 

T)  “  0:  g  =  0  ,  g  *  0  ,  V  *  0  (62) 

1  *  8  «  g»  V  «  1  ,  h  *  1  . 

(63) 


A  numerical  example,  which  is  determined  by  the  Reynolds  number 
R  *  10,  is  displayed  in  figure  2.  The  solution  has  been  obtained 
by  the  Runge-Kutta  method  which  started  the  integration  with  an 
assumed  set  of  initial  data  that  had  to  be  improved  successively 
until  the  accuracy  conditions  (B)  and  (c)  (see  section  3)  were 
sufficiently  met. 

A  more  efficient  iteration  procedure,  which  improves  Successively 
the  crude  solution  of  the  system  (51)  and  (52)  by  integrating  the 
linear  equations 


8n  +  Tfe„  -  f  (gU  +  V2  -  h^) 
V  +  T]V  «  0 

8,(8r.  -  Hg,)  -  -i  hn  , 


(64) 

(65) 


(66) 


will  be  described  in  another  paper  in  preparation.  Numerical 
results  will  be  presented  and  discussed  for  various  Reynolds 
numbers  R. 

All  numerical  calculations  reveal  a  very  strong  dependence  of  the 
secondary  flow  (u,  w)  upon  the  Reynolds  number  R.  Hence,  it  is  useful 
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to  introduce  a  "boundary  layer  of  the  secondary  flow"  the  thickness  of 
which  Will  depend  upon  the  Reynolds  number  R.  This  boundary  layer  may 
be  determined  by  the  parameter  value  Q  such  that  at 

C  =  Cs  :  W(0  m%-4  .  (67) 

The  limiting  line  z  »  6g(r)  of  the  boundary  layer  of  the  secondary 
flow  is  then  defined  by 

z  =  6s(r>  s  Cs6(r)  , 

where  z  =  6(r)  is  the  limiting  line  of  the  boundary  layer  of  the 
entire  flow  defined  by  the  similar  condition  (47)  (see  figure  1  and  2). 

After  simple  numerical  calculations  one  arrives  at  the  following 
two  significant  properties  of  vortex  flows. 

Phenomenon  5:  The  boundary  layer  thickness  of  the  secondary  flow 
decreases  as  the  Reynolds  number  R  increases. 

Phenomenon  6:  At  the  solid  surface  the  radial  shear  stress  increases 
much  faster  with  the  Reynolds  number  R  than  the  tangential  shear  stress. 

The  phenomena  1  through  6,  which  appear  to  be  new  in  the  theory  of 
real  flows,  are  physically  plausible.  This  has  been  pointed  out  for 
the  phenomena  1  and  2  in  the  previous  sections.  The  phenomena  3 
through  6  are  physically  also  feasible,  as  an  increasing  tangential 
velocity  tends  to  decrease  and  to  increase  the  boundary  layer  thickness 
at  the  same  time,  indeed,  when  the  particles  of  the  fluid  are  forced 
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to  remain  closer  to  their  undisturbed  circular  motion  they  are  simul¬ 
taneously  forced  to  remain  longer  under  the  influence  of  the  friction 
forces  at  the  surface.  Thus ,  an  increasing  tangential  velocity  can 
neither  increase  nor  decrease  the  thickness  of  the  friction  layer  at  the 
surface.  However ,  these  arguments  do  not  apply  to  the  secondary  flow. 
Hence  the  boundary  layer  of  the  secondary  flow  must  behave  in  the  usual 
manner ,  i.e. ,  its  thickness  must  decrease  as  the  Reynolds  number 
increases . 

The  properties  of  vortex  flows  over  a  flat  surface  found  above 
may  be  compared  with  phenomena  observed  outside  the  cores  of  hurricanes. 
Such  a  comparison  is  generally  feasible  ,  as  the  driving  core  of  a 
hurricane  may  roughly  be  replaced  by  a  rotating  rod  that  produces  a 
similar  vortex  flow.  Of  course  exterior  disturbances  of  the  vortex 
flow  other  than  those  caused  by  the  surface  of  the  earth  must  be 
neglected.  The  observations  have  fully  confirmed  the  source  and 
sink  character  of  the  boundary  layer  at  the  surface  for  the  secondary 
axial  flow.  The  logarithmic  increase  of  the  axial  velocity  toward 
the  vortex  axis  exp lains  the  rapidly  increasing  rainfall  toward  the 
core  of  a  hurricane  (see  [10,  p.  130  ) .  In  this  connection  it 
appears  feasible  to  identify  the  radius  of  the  ra in  fa 1 1  area  roughly 
with  the  characteristic  flow  parameter  r0 ,  which  measures  the  radial 
extent  of  the  logarithmic  source  around  the  axis.  The  radial  shear 
stress  at  the  surface,  which  is  large  compared  with  the  tangential 
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shear  stress ,  explains  the  very  strong  radial  ocean  waves  produced  by 
hurricanes  (see  [6,  p.  298])*  Finally,  it  may  be  mentioned  that  the 
temperature  observed  in  hurricanes  (see  .6 ,  p.  319 _)  exhibits  a  similar 
nonmonotonic  behavior  as  the  axial  velocity  found  above.  This  should 
be  expected,  because  the  energy  equation,  which  governs  the  tempera¬ 
ture  of  the  fluid,  has  the  same  second  order  terms  as  the  third 
Navier-Stok.,3  equation  for  the  ax' a’  velocity 
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FIGURE  2  :  The  dimensionless  velocity  (U,  V,  W)  and  the  pressure 
2P  vs.  the  dimensionless  variable  ?.  for  the  Reynolds 
number  10. 
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